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ABSTRACT 

We report deep near-infrared and optical observations of the X-ray point 
source in the Cassiopeia A supernova remnant, CXO J232327. 9+584842. We 
have identified a J = 21.4±0.3 mag and Kg = 20.5 ± 0.3 mag source within 
the l-cr error circle, but we believe this source is a foreground Pop II star with 
Tgfj = 2600-2800 K at a distance of ~ 2 kpc, which could not be the X-ray point 
source. We do not detect any sources in this direction at the distance of Cas A, 
and therefore place 3-0" limits of i? > 25 mag, F675W > 27.3 mag, J > 22.5 mag 
and Kg > 21.2 mag (and roughly if > 20 mag) on emission from the X-ray point 
source, corresponding to Mji > 8.2 mag, Mpeysw ^ 10-7 mag, Mj > 8.5 mag, 
Mh > 6.5 mag, and Mk^ ^ 8.0 mag, assuming a distance of 3.4 kpc and an 
extinction Ay = 5 mag. 

Subject headings: infrared: stars — stars: late-type — supernova: individual 
(Cassiopeia A) 

1. Introduction 



Cassiopeia A (Cas A) is the youngest Galactic supernova remnant (SNR) with an age 
of ~ 320 yr, as according to Ashworth (1980) it is associated with the explosion observed by 
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Flamsteed (1725) in 1680. Hughes et al. (2000) have found that the elemental abundances 
in Cas A are consistent with those expected from the remnant of a massive star, possibly a 
Wolf-Rayet star (Fcsen et al. 1987), and therefore Cas A is considered to have been a Type II 
supernova. One therefore expects a compact central remnant, such a neutron star or black 
hole, based on the initial mass function of Type 11 supernovae (e.g. dc Bonder & Vanbcvcrcn 

1998) . From the first-light images of the Chandra X-ray Observatory (CXO), Tananbaum 
(1999) reported detection of a compact source located at the apparent center of Cas A. The 
detection of this source, CXO J232327.9+584842 (hereafter the X-ray point source or XPS), 
was later confirmed in archival ROSAT (Aschenbach 1999) and Einstein (Pavlov & Zavlin 

1999) data. 

The XPS is located within 5" of the expansion center of Cas A (van den Bergh & 
Kamper 1983), and given the space density of AGN the chance of finding one within this 
distance of the center is quite small. We convert the count rates from Chakrabarty et al. 

(2001) to the 0.5-2.4 kcV band, and get an absorbed flux of ~ 4 x lO"-*-^ ergs s~^ cm~^. 
Comparing this with the AGN logiV-logS" relation from Georgantopoulos et al. (1996), we 
would expect -^0.4 AGN deg"^ or ~ 2 x 10"® AGN of this flux at the center of Cas A. 
It is thus extremely improbable that the XPS is an AGN, a fact further confirmed by its 
relatively steep spectrum (Chakrabarty et al. 2001). 

Therefore, it is generally believed that the XPS is associated with the remnant of the 
Cas A progenitor (Chakrabarty ct al. 2001). The X-ray spectrum of the XPS, as determined 
by Pavlov et al. (2000) and Chakrabarty ct al. (2001), can be fitted by a power-law with a 
photon index ~ 3. Other acceptable fits include thermal bremsstrahlung [kT^ 1.7 keV), 
blackbody {kT°° f« 0.5 keV, R°° ^ 0.5 km), or neutron star atmospheres {kT°° 0.4 keV, 
R°° ^ 0.8 km for the model of Heyl & Hernquist 1998; kT"^ ^ 0.27 keV, 2 km for the 

model of Zavhn et al. 1996). 

The nature of the XPS is unclear. However, we have an idea as to what it is not. The 
spectral index of the XPS is significantly steeper than those typical for young X-ray pulsars, 
its luminosity is > 10^ times less than those of young X-ray pulsars, and there is no evidence 
for a synchrotron nebula (McLaughhn et al. 2001). The spectrum is similar to that of an 
anomalous X-ray pulsar (AXP; see Mereghetti 2000), but the X-ray luminosity is at least 
several (if not 10-100) times fainter than that typical for AXPs. The XPS is cooler but 
much more luminous than isolated neutron stars (Motch 2000). 

Furthermore, there have not been any detections of optical (van den Bergh & Pritchet 
1986; Ryan et al. 2001) or radio (McLaughlin ct al. 2001, and references therein) emission 
from the XPS, nor have X-ray pulsations been detected (Chakrabarty ct al. 2001), though 
the current limits are not very constraining. Therefore, the XPS is almost certainly not a 
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young pulsar similar to the Crab. Theories as to its identity range from a cooling neutron 
star emitting from polar caps to an accreting black hole (Umeda et al. 2000; Pavlov et al. 
2000). 

Prom measures of hne ratios in the Cas A remnant, Searle (1971) finds the extinction 
to be Ay — 4.3 mag in the direction of one filament. Later radio studies found significant 
variations of Ay on scales of ~ 1', and overall values ranging from 4-5 mag for the north and 
northeastern rim and > 5-6 mag for the rest of the SNR (Troland ct al. 1985). Similarly, 
Hurford & Fcscn (1996) find extinction values of 4.6-5.4 mag across the northern portion 
(assuming Ry = 3.1). We will therefore adopt a middle value of Ay ~ 5 mag. We assume 
that Cas A and the XPS arc at a distance of 3.4^q;J kpc (Reed et al. 1995), which we 
parameterize as D — 3.4^3 4 kpc. 

In this letter we report on optical/near-IR searches for a counterpart to the XPS. Wc 
believe that given the unknown nature of XPS, searches at all wavelengths arc warranted and 
even upper limits may constrain the nature of this enigmatic source. The paper is organized 
as follows: in Section 2 we detail our observations and reduction techniques. Section 3 
contains a description of the results, while Section 4 presents an analysis of these results. 
Finally, a discussion and conclusions are in Section 5. 

2. Observations 

2.1. Cas- A Central Point Source Position 

The SNR Cas A was observed several times with the CXO. After the initial detection 
in the first-light images (Tananbaum 1999), a long HRC-I observation was obtained on 1999 
December 20, and a third observation with the HRC-S in imaging mode was taken on 2000 
October 5. A discussion of the results from this observation is in preparation (Murray et 
al. 2001); here we provide only the source location information. Table 1 gives the point 
source locations and estimated uncertainties (including estimates of systematic errors). We 
estimate that the overall positional uncertainty for all of these observations is I'.'O (1-cr). 

2.2. Optical and Neeir-IR Observations 

The observations were carried out primarily with the Near Infrared Camera (NIRC; 
Matthews & Soifer 1994) mounted on the 10-m Keck I telescope, augmented with archival 
HST/WFPC2 images. We also took auxiliary optical and infrared calibration images with 
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the COSMIC imager on the Palomar 5-m telescope (P200), the PFIRCAM (Jarrett et al. 
1994) infrared imager on the P200, and the P60CCD optical imager on the Palomar 1.5-m 
telescope (P60). A summary of the instruments, filters, exposures, and conditions is listed 
in Table 2. 

The optical data were reduced with the standard IRAF ccdred package. The images 
were bias subtracted, flat-fielded, registered, and co-added. The infrared data were reduced 
with custom IRAF software. The images were dark subtracted, fiat-fielded, and corrected for 
bad pixels and cosmic rays. We then made object masks, which were used in a second round 
of flat-fielding to remove holes from the flats. The data were then registered, shifted, and 
co-added. The HST images were processed using the standard drizzling procedure (Fruchter 
& Mutchler 1998). 

The data from the P60CCD were used as the astrometric reference. We matched 36 
non-saturated stars to those from the USNO-A2.0 catalog (Monet 1998). Using the task 
ccmap we computed a transformation solution, giving 0'/2 residuals (all astrometric residuals 
are l-cr for each coordinate unless otherwise indicated). We then used this solution to fit 
stars on the COSMIC images. Using 37 stars, we again obtained (X'2 residuals. 

We then used 15 stars on the COSMIC images to fit the HST image, getting (X'O? residu- 
als. This solution was then used for the infrared images, fitting 10 stars on the NIRC images 
with (y.'OS residuals. This gives (y/4 position uncertainties relative to the ICRS, assuming the 
uncertainties intrinsic to the USNO-A2.0 are (y/3 (for each axis; Monet 1998) We then trans- 
fered this solution to the PFIRCAM images (0'.'3 residuals), but as this is only a photometric 
reference the absolute position is not important. 

For the optical photometry, we used V, R, and / observations of the standard fields^ 
Landolt 110, NGC 7790, and PC 1657 (Landolt 1992; Stetson 2000) carried out with the 
P60CCD. We fit the observations over the whole night using airmass corrections and first- 
order color terms, and measured the R zero-point magnitude. We then examined 25 stars 
on the Cas A images common to both the P60CCD and COSMIC images, and from this 
determined the zero-point for the photometric night. From these data we also determined 
the limiting magnitude tobei?~25 mag. 

For the infrared photometry, we used 3 observations of the faint UKIRT standard stars 
FS 29 and FS 31 (Casah & Hawarden 1992) taken with the PFIRCAM. These observations 
were used to determine J, H , and Kg zero-points (we assumed the magnitudes were the 
same as the K magnitudes, as the correction is typically < 0.01 mag: much smaller than 



^http://cadcwww.dao.nrc.ca/cadcbin/wdb/astrocat/stetson/query/ 
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our uncertainties; Persson et al. 1998). Prom these images we then found 5 stars common 
to the PPIRCAM and NIRC images, and determined zero-point magnitudes for NIRC. 



3. Results 

We searched for a counterpart to the X-ray point source, at position q;(J2000) = 
23'^23'"27!857, 5( J2000) = +58°48'42'.'77, with I'/O uncertainty (Table 1). See Pigure 1 for the 
separate optical/IR images. There was no source on COSMIC images, giving R > 25 mag 
(3-0" limit) for any possible counterpart (this agrees with the previous limit of i? > 24.8 mag 
and I > 23.5 mag; van den Bergh & Pritchet 1986). On the NIRC, PPIRCAM, and HST 
images there was a source 1'.'7 away from the X-ray position, at «(J2000) = 23'^23™27!78, 
5(J2000) = -|-58°48'41'.'2 (±0'.'4 in each coordinate). Given the astrometric uncertainties, the 
overall position uncertainty is I'.'l in each axis, so this source is 1.5-a away from the nominal 
position. We label this source A, and consider it as a potential candidate counterpart or 
companion to the X-ray source. The magnitudes of source A are P675W = 26.7 ± 0.2 mag, 
J = 21.4 ± 0.3 mag, H ^ 20.5 ± 0.8 mag, and = 20.5 ± 0.3 mag. There are no other 
sources within the 2'/3 radius 90% confidence circle. 



4. Analysis 

Using the reddening and zero-point cahbration data from Bessell et al. (1998), we plot 
the spectral energy distribution (SED) for source A in Pigure 2. This incorporates both the 
detections and limits. 

To determine if source A could be a star, we compared model stellar colors from Bessell 
et al. (1998) with our data. We fitted for three parameters: the visual extinction Ay, the 
distance in kpc Dkpo ^-nd the stellar model (which includes the effective temperature Teflf, 
the surface gravity g, and the metaUicity [Pe/H]). We assumed that the star would be a 

zero-age main-sequence star such that log{R/ Rq) = O.71og(M/M0) —0.1 (Habets & Heintze 
1981), and used the bolometric corrections and reddening from Bessell et al. (1998) to find 
the expected magnitudes. To account for the upper limits in our fitting, we minimized a 
modified statistic, such that 



/rn -m- ^ \ ^ non-Detect f Q if mj,mod > m; 

^ - 2^ [ a, J+Z^l (^ ra,,^o,-m, ^^ OthcrwiSC 



Here, i runs over the different filters, mj is the observed magnitude or limit for that filter. 
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?7ii,mod is the model magnitude, and (Tj is the uncertainty. The model uses standard Vega- 
based magnitudes, where the HST data do not. Therefore we converted the HST magnitude 
to the Vega-based system using R — F675W = —1.05 mag, appropriate for sources of this 
color^. We do not incorporate model uncertainties into this statistic. Minimizing this seeks 
the best model that comes close in magnitude to the detections while remaining fainter than 
the non-detections. A full-fledged Bayesian analysis (e.g. Gregory & Loredo 1992; Cordes & 
Chernoff 1997) would be more accurate, but we only wish to demonstrate the plausibility of 
model fits, not assign specific probabilities to different models. 

Given the number of variables, this fit is somewhat unconstrained. We restrict the 
extinction and distance to reasonable values (0.5 < Ay < 8 mag, 0.5 < -Dkpc ^ 5), and 
fit for log(5i) = 5.0 (appropriate for late M stars; Habets & Heintze 1981). In addition, we 
require that Ay roughly scale with excluding models that are very distant but have almost 
no extinction. We find that our detections and hmits are entirely consistent with a cool (M6- 
8), Pop II, main-sequence star, which is between the Earth and Gas A. Good fits are obtained 
for stars with %s = 2600-2800 K, [Fe/H]= -2.0, Dkpc = 1.8-2.0 kpc, and Ay = 3.1-3.2 mag 
(see Figure 2 for examples). We do not give the value or formal confidence regions, as 
the ill Equation 1 is somewhat contrived and the models for stars this cool are not well 
determined (Bessell et al. 1998), but Figure 2 demonstrates the plausibility of the fits. That 
there is a star within 1'.'7 of the XPS is quite behevable: the theoretical star-count models of 
Nakajima et al. (2000) give 1.5 x 10^ stars deg~^ of the appropriate colors with J < 22.5 mag, 
leading to a false coincidence rate of 1.0. The best-fit star has R — 0.2Rq, M — O.IMq, and 
L — O.OO4L0. Shghtly deeper I band observations should be able to verify the classification 
of source A. 

As one might expect, there is a significant anti-correlation in the fits between values of 
£)kpc and Ay, with ±0.25 kpc and ±0.5 mag variations giving reasonable fits, but the fits 
for the range of hkely extinctions for Gas A (4-6 mag) at 3.4 kpc are definitely poor. 

Assuming that source A is a late-type star, we examine whether it could be associated 
with the XPS, implying that both are in the foreground and that the XPS is not associated 
with Gas A. Prom Katsova & Gherepashchuk (2000), we see that for a star with B — V > 1.8 
(from the model for source A), the X-ray luminosity is Lx,star(0. 1-2.4 keV) < 10^^ ergs s~^ 
(also James et al. 2000; Marino et al. 2000), giving unabsorbed (denoted by superscript U) 
X-ray-to- infrared flux ratios of fx/{^jfu,j) = 6 x 10"^ and f^/ {i^Ksfu,Ks) = 2 x 10"^ for such 
a star. Gonverting the flux of the XPS to the ROSAT passband (using W3PIMMS3), it has 



^http:/ /www. stsci.edu/instruments/wfpc2/Wfpc2_phot/wfpc2_cookbook. html 
^http://heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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ratios of fx/{^jfv,j) — 34, fx/{^Ksfv,Kj — which are drastically different. In addition, 
Pavlov et al. (2000) and Chakrabarty et al. (2001) did not observe any variability from the 
XPS, unlike late-type stars that can vary by factors of ~ 10^ on small time scales (Marino 
et al. 2000). Source A therefore could not emit the X-rays observed from the XPS. 

We conclude that the XPS was not detected, and add F675W > 27.3 mag, J > 22.5 mag 
and Kg > 21.2 mag (S-cr), along with a rough limit of if > 20 mag, to the previously 
mentioned limits. 



5. Discussion & Conclusions 

Based on a synthesis of CXO, ROSAT, and Einstein data, Pavlov et al. (2000) fit the 
X-ray spectrum of the XPS. The absorbed flux is 8.2 x 10~^^ ergs cm~^ s~^ in the 0.3- 
6.0 keV range (Pavlov et al. 2000). Power-law and pure blackbody models give good fits to 
the absorption-corrected data, and are plotted as representative X-ray spectra in Figure 3. 
These results are similar to those from Chakrabarty et al. (2001). Pavlov et al. (2000) prefer 
the results of a H/He polar-cap model with a cooler Fe surface, but all we wish to illustrate is 
that blackbody models are consistent with the optical limits, while power-law models require 
a break between the X-ray and optical bands. 

In Figure 3 we also plot the expected optical magnitudes of representative X-ray sources 
(an AXP and a tight X-ray binary) for comparison. These magnitudes are derived by taking 
the X-ray-to-optical flux ratios for these objects and scaling them to the X-ray flux of the 
XPS. We can hkely reject sources hke 4U 1626—67 (Chakrabarty 1998) from consideration, 
but the extrapolation of the AXP 4U 0142-1-61 (HuUeman et al. 2000) is consistent with the 
current limits. 

Giving the presumed distance and reddening, our limits translate to Mr > 8.2 mag, 
Mf675W > 10.7 mag, Mj > 8.5 mag, Mh > 6.5 mag, and Mk, > 8.0 mag. We find the 
observed X-ray-to-infrared fiux ratios to be /x/(i^F675w/i.,F675w) Z 2872, fx/{i^jU,j) > 212, 
f^l {^Ksfv,Ks) ^ 280 (the X-ray flux is in the 0.3-6.0 keV band). If we correct for interstellar 
absorption, we find unabsorbed ratios of /x/('^F675w/"f675w) ^ 231, fx/{^jfu,j) ^ 166 and 
fx/i^KsfuKs) ^ ^67, using the X-ray flux from Chakrabarty et al. (2001). These flux ratios, 
larger than those inferred previously, tighten constraints on the identity of the XPS (e.g. 
Umeda et al. 2000; Pavlov et al. 2000). 
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Table 1. Cas A X-ray Observation Summary 



OBSID 


Date 


CXO 
Instrument 


Exposure 
(ks) 


RA^ 
(J2000) 


Dec* 
(J2000) 


214 


1999 Aug 20 


ACTS S3 


6 


23'^23"^27!94 


+58°48'42'/4 


1505 


1999 Dec 20 


HRC-I 


50 


23''23"'27^88 


+58°48'42'/l 


1857 


2000 Oct 05 


HRC-S 


50 


23''23"'27^75 


+58°48'43'/8 


Average 
Uncertainty'' 








23''23"'27!857 
ffi097 


+58°48'42'/77 
0'/91 



*Thc individual source positions were calculated as centroids of the event distributions 
taken within a I'.'O radius circle about the location and iterated until the centroid location 
shifted by less than O'.'l. 

^Uncertainties are 1-a. 



Table 2. Cas A Optical/Near-IR Observation Summary 



Date 


Telescope / 
Instrument 


Observer 


Band 


Exposure 
(s) 


Conditions 


2000 Jan 22 


HST/WFPC2 


R. Fesen 


F675W 


4000 




2000 Jun 27 


Keck I/NIRC 


S. Kulkarni 


J 


1600 


slight cirrus 








Ks 


2364 




2000 Jul 04 


P200/COSMIC 


P. Mao 


R 


1010 


photometric 


2000 Jul 05 


P200/COSMIC 


P. Mao 


R 


1000 


high cirrus 


2000 Jul 24 


P60/P60CCD 


D. Kaplan 


R 


150 


photometric 








I 


150 




2000 Sep 06 


P200/PFIRCAM 


D. Kaplan / 


J 


1680 


photometric 






J. Cordes 


H 


1680 










Ks 


1120 
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Fig. 1. — Images of the region around the XPS. They are: i?-band (COSMIC; upper left); 
J-band (NIRC; upper right); //-band (PFIRCAM; lower left); X.-band (NIRC; lower right). 
A 2'.'3 radius circle (90% confidence) is drawn around the position of the XPS, and candidate 
source A is indicated. North is up, and East is to the left. The images are f« 15" on each 
side. 
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V (Hz) V (Hz) 

Fig. 2. — Spectral energy distribution for source A, a foreground star from the Cas A X- 
ray error circle, and best-fit data for two different stellar models: T^s = 2600 K (left); 
Teff = 2800 K (right). The open symbols are the observed data, the filled circles and limits 
those corrected for reddening, and the diamonds the model values. The model parameters, 
from Bessell et al. (1998), are hsted on the figures. 
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Fig. 3. — Spectral energy distribution for the XPS. This incorporates optical hmits (this 
work) and X-ray data (Pavlov et al. 2000). The open triangles are the measured values, 
while the filled triangles are those corrected with Ay = 5.0 mag. The open circle is the 
measured CXO flux, while the thick lines are model spectra corrected for absorption: power- 
law (PL; dashed) and blackbody (BB; dash-dotted). The thin dashed lines are derived from 
the unabsorbed X-ray-to- i?-band ratios of the AXP 4U 0142-^61 (HuUeman et al. 2000, 
with Ay — 4:A mag) and the very close X-ray binary 4U 1626—67 (Chakrabarty 1998, with 
Av — 0.2 mag), assuming a 0.5-10 keV luminosity of 10^^ ergs s~^ for the XPS. We do not 
plot the more comphcated atmosphere models from Pavlov et al. (2000) or Chakrabarty 
et al. (2001). 



